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We examine the effect of superluminal signal propagation through a birefringent crys- 
tal, where the effect is not due to absorption or reflection, but to the filtration of a special 
polarization component. Wc first examine the effect by a stationary phase analysis, with 
^ (— I results consistent with those of an earlier analysis of the system. We supplement this 

analysis by considering the transit of a gaussian wave and find bounds for the validity 
j of the stationary phase result. The propagation of the gaussian wave is illustrated by 
figures. 

:3 

lH]^ Introduction 

For wave propagation in a medium, group velocity is often regarded as being identical to signal 
J — velocity, i.e., the velocity at which energy, and thereby information, is transferred. However, 

lO it has been known for a long time that group velocity, as usually defined, can exceed the 

speed of light without breaking causality [1, 2, 3]. In more recent years such superluminal 
propagation has been in the focus of interest, in particular since a series of different types 
I of experiments now have been performed to demonstrate the effect. These experiments have 

been done in media with anomalous dispersion and under conditions where superluminal 
tunneling times of single photons have been observed [4, 5, 6, 7]. 
I> Anomalous dispersion, and thereby superluminal group velocity, is usually obtained in 

systems with absorption or reflection. The relation between damping and group velocity 
^ is seen most directly in linear systems, where the Kronig-Kramers relations link the real 

part of the propagator to the imaginary part. The peculiar fact is that the KK relations, 
which on one hand reflects the condition of causality, on the other hand shows explicitly how 
superluminality arises in systems with anomalous dispersion. 

Thus, even if group velocities in these systems may exceed the speed of light in vacuum, 
there is no conflict with causality. The apparent paradox, superluminality without breaking 
of causality, can be viewed as a consequence of the definition of group velocity, Vg = duj/dk. 
In systems with anomalous dispersion Vg typically varies rapidly with frequency lu in regions 
with large Vg, and to interpret this as the velocity of propagation therefore is meaningful only 
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for signals that are sharply defined in frequency and therefore slowly varying in space and 
time. A signal with a sharper space profile will typically be rapidly deformed and damped 
so that superluminal propagation speed cannot be obtained over long distances. In fact for 
signals with a well defined front, the front velocity can never exceed the speed of light, and 
this fact is also the essence of the KK-relations. 

Even if anomalous dispersion cannot occur in a transparent, passive medium, it is pos- 
sible in an active medium with gain-assisted wave propagation. This has been shown both 
theoretically and experimentally [8, 7]. In such a case the signal is not damped, but the 
theoretical group velocity will also here only in an approximate sense represent the real speed 
of propagation for a signal with finite space extension. Thus, as discussed in [8, 9] the signal 
propagation can be viewed as a pulse reshaping process. In this process the signal peak may 
move faster than light, but only for a limited time, since the signal front cannot be overtaken. 

Some time ago D.R. Solli et.al. showed that even in a passive transparent medium anoma- 
lous dispersion can be achieved by filtration of a component of a signal with normal dispersion 
[10, 11]. They demonstrated this, theoretically and experimentally, for a strongly birefringent 
photonic crystal with electromagnetic waves of centimeter wave length. With an incoming 
wave that decomposes equally in a fast and a slow component and with an outgoing wave 
that is filtered in the same polarization direction as the incoming signal a dispersion relation 
was found that under half waveplate conditions could be interpreted as superluminal transit 
speed through the crystal. This result is interesting, since on one hand the anomalous dis- 
persion of the filtered signal clearly is consistent with the interpretation of a superluminal 
group velocity and on the other hand it is clear that in this case a faster-than-light effect can 
only be apparent. Thus, the signal is created as a linear superposition of two waves, each 
of which are transmitted with subluminal velocity through the crystal. (In the context of 
quantum "weak value" measurements, the effect of the polarization filter can also be seen as 
performing a "postselection" of the photon state, as previously discussed in Ref. [12].) 

The results of ref. [10] is the motivation for the present work, where we examine how the 
derived group velocity is related to the transit of wave packets through the crystal. We first 
reproduce the earlier results by use of a stationary phase analysis and focus then on how a 
signal of gaussian shape is moving through the crystal. The peak position of the outgoing 
signal is found to correspond to superluminal transit only when the wave packet is much 
wider than the crystal, and even in this case the filtered signal lies well within the envelopes 
of the two subluminal components, corresponding to the fast and slow wave components in 
the crystal. 

Plane wave propagation and peak velocity 

The physical system we consider is sketched in Fig. 1. A polarized electromagnetic wave is 
sent through a birefringent crystal, which is oriented with the optical axis (z-axis) orthogonal 
to the direction of propagation (x-axis). The boundary planes of the crystal, where the 
wave enters and leaves the crystal, are both assumed to be orthogonal to the direction of 
propagation. The outgoing wave passes a polarization filter so that the registered signal 
corresponds to linear polarization at a rotation angle P in the y,z-plane. We consider here 
how the amplitude of the filtered signal depends on the angle /3 and the wave number k of 
the incoming wave. 

The crystal is assumed to be linear and non- absorbing, with a dielectric tensor that is 
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Figure 1: Schematic illustration of the experimental situation considered in the text. A plane 
polarized electromagnetic wave is sent through a birefringent crystal C of width d. After 
leaving the crystal the wave passes a polarization filter F with a polarization angle (3 that 
can be varied, and the amplitude of the filtered signal is registered in the detector D. 

frequency independent. Due to the difference between index of refraction ni outside and 712 
inside the crystal, the reflection coefficient [13] 



is slightly different from 0. In fact, since the index of refraction is different for the ordinary 
wave (n2 = no), with polarization along the y-axis, and the extraordinary wave (n2 = ne), 
with polarization along the z axis, the reflection coefficient for these two directions of polar- 
ization are not identical. However, for both directions we will consider the deviation from 
to be sufficiently small so that the crystal can be treated as non-reflecting. (In Ref.[ll] 
the average index of refraction is given as about 1.25 corresponding to a reflection coefficient 
R ~ 0.01.) We shall in the following assume that An = — Uq to be positive so that the 
ordinary wave is the fast and the extraordinary wave is the slow component in the birefringent 
material. 

The incoming wave is assumed to be linearly polarized, with direction of polarization ro- 
tated 45° relative to the y-axis. This means that the wave, inside the crystal, has components 
of ordinary and extraordinary waves of equal amplitudes. Due to the different propagation 
velocity of the two components in the crystal, the polarization will change to elliptical polar- 
ization during the propagation. For a distance corresponding to half-waveplate condition the 
polarization is again linear, with polarization direction orthogonal to the polarization of the 
incoming wave. In the following we assume the length d of the crystal in the direction of propa- 
gation to be close to that determined by the half waveplate condition Ak d = (2A^-|-l)7r, N ^"L, 
with AA: = — ko as the difference between the wave number of the extraordinary and or- 
dinary wave. This means that the wave which leaves the crystal will be close to linearly 
polarized in the direction orthogonal to that of the incoming wave. The filtered signal for 
a polarization direction (3 close to that of the incoming signal therefore picks out only a 
small component of the outgoing wave, but this is the interesting one for the question of fast 
propagation through the crystal. 

To analyze the situation, we start with the expression for the electric component of a 
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monochromatic plane wave propagating through the crystal 



E(x.t) . ^ 



^i{k„X-Ujt)^ _^ ^i{keX-0jt)2 0<X<d (2) 

^i{kx-ujt) ^^i(ko-k)d^ _j_ ^i{ke-k)d£^ d < X , 



where y and z are the unit vectors in the y and z directions. The three intervals on the x-axis 
correspond to the regions where the wave is located before entering the crystal (x < 0), when 
inside the crystal < x < d) and after leaving the crystal {d < x). The polarization filter 
at this point has not been introduced. The wave numbers in the three regions are related to 
the frequency by 

K = , Ko ~ '^O) f^e ~ '^e ) 
C C C 

with Uo and Ug as the indexes of refraction of the ordinary and extraordinary waves, respec- 
tively. 

The filtered signal, with polarization at angle (3 relative to the y-axis is then described by 
the electric field 

Ef) = -E-ef3= ^e^('=^-'^*)e'('=°-*^)'^(cos /3 + sin /3 e^^'='^) , (4) 
with A A; = kg — ko- We focus the attention on the complex, relative amplitude 

z= -L(cos/3 + sin/3e^^'='^) = |z|e^>^, (5) 

with absolute value and complex phase 

\z\ = ^[1 + sin 2/3 cos{Akd)]/2 

/ sm(Akd) \ (6) 

Y = arctan , . , . 

\cos{Akd) + cot P J 

The amplitude \z\ vanishes at isolated points in the {k, /3)-plane, determined by the conditions 

a) /3 = 7r/4 (57r/4), Afc d = tt (mod 27r) 

b) 13 = Sir /A (77r/4) , Akd = (mod 2tt) , (7) 

with a) corresponding to the half waveplate condition. The phase x rotates by 27r around 
each point of zero amplitude, and the rapid variation close to this point is the reason for 
the large value of the group velocity near the half wave plate frequency. In Fig. 2 a graphic 
representation of the amplitude z is shown as a function of k and /3 in the form of contour 
plots for \z\ and x- 

The expression (6) for the phase x of the filtered plane wave contains information not 
only about plane wave propagation in the crystal, but also, in an approximate sense, about 
propagation of wave packets, which in k-space are strongly localized around the wave number 
of the plane wave. Thus, a standard stationary phase analysis is based on the assumption 
that the peak position of the wave in x-space is determined by the condition that the complex 
phase of the plane wave is stationary with respect to variations in k. 
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Figure 2: Amplitude and phase of the filtered monochromatic plane wave as functions of the 
polarization angle /? and the wave number k. The (red) curves that circulate the singular 
points are contour lines of constant amplitude and the (blue) lines that interconnect these 
points are contour lines of constant phase. Three singular points with vanishing amplitude 
are shown, where the lower one corresponds to half waveplate condition for the transmitted 
wave. The numbers along the horizontal axis gives A/cc? measured from the closest point with 
an integer number of wave lengths {2tiN). 



For the filtered wave the time dependent position of the peak is then given by 



{kx -Ljjt + {ko - k)d + x) = 



X 



ct + d — Uod — c 



dio 



(8) 



Clearly the stationary phase argument identifies the group velocity, defined in the usual 
way V group — ^ the velocity of the peak of the signal. However, when we consider the 
propagation of the wave after leaving the crystal there is no change relative to the velocity 
of light in vacuum. More interestingly, the expression gives the time delay (or time advance) 
relative to a wave that moves freely (without the presence of the crystal) [14] . The shift in 
time of the outgoing, filtered signal relative to a freely propagating signal is 



dx ^ {no - 
du c 



l)d 



(9) 



and with x given by (6) we find 

1 + cos(Afcd) cot P 



d 



1 + cot^ /5 + 2 cos(A/cd) cot /3 c 
cos 2/3 An 



n„ 



+ 



{uo - l)d 



d 



n 



1 



1 



2 1 + cos(A/c(i)sin2/3 



(10) 



where in the last expression we have introduced n = {Uf, + n^jl and An = Tie 

Eq.(9) shows that the time shift r becomes negative when ^ is sufficiently large and 
negative. 



duj 



< 



{no - l)d 



(11) 
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Figure 3: Lines of constant time shift for the transmitted signal in the /?, fe-plane. The 
lower (red) curves correspond to advanced signals and the upper (blue) curves correspond to 
retarded signals. Close to the point of vanishing amplitude the time shift increases without 
any bound. The region of superluminal transit, for parameter values given in the text, is 
indicated by shaded red. 

This condition, which one may interpret as the condition for superluminal transit through the 
crystal, can be satisfied for parameter values close to the points of vanishing amplitude, z = 0, 
where the phase angle is rapidly varying. To be more specific, we introduce new variables 93 
and e. 



and assume these to be sufficiently close to the point ip = e = so that an expansion to first 
order in these variables is sufficient. 



This gives for the derivative of the phase X: 

dx 2e 1 _ 2e I d 

duj ~ Me2 + (^2 + 2U0J ~ Me2 + v?2 + 2^c ' 

and the corresponding expression for the time shift of the peak, valid for small 99 and e. 




(12) 




(13) 




(15) 



The expression shows that for negative e and for sufficiently small and ip^^ the time shift r 
may become negative. In fact, close to the singular point of the phase x the time shift may 
become arbitrarily large and negative. Thus, the transit time may not only correspond to a 



superluminal propagation, but to a situation where the peak of the outgoing, filtered signal 
appears before the peak of the incoming signal has entered the crystal. Even if this situation 
may appear paradoxical, there is in reality no conflict with causality. The explanation is 
simply that the outgoing signal is being produced by the advanced tail of a sufficiently broad 
incoming wave. 

The condition for super luminality (11) for small e and reduces to the condition 

(e + a)^ An , ^ 

a2 Aa? 4(n - 1) 

In Fig. 3 the contour lines of constant time shift r are shown around the point e = </? = 0. 

The time shift is positive for /3 > -7r/4 (e > 0) and negative for /3 < 7r/4 (e < 0), with 
largest absolute values of r close to the singular point e = = 0. In the figure the region of 
superluminality is indicated by the shaded area for parameter value a = 0.15. 

The results discussed so far are in full agreement with those discussed in Ref.[10]. However, 
these results are based on the use of the stationary phase argument, which can only be trusted 
for a signal that is narrow in k-space and therefore wide in x-space. In the next section we 
shall therefore supplement this discussion by examining the propagation through the crystal 
of gaussian shaped electromagnetic waves of different widths. 



Propagation of Gaussian shaped wave packets 



We introduce a wave packet with gaussian envelope in the direction of propagation by the 
following expression 



(17) 



E(a;,t) = j c/a;e-j('^-'^°)'E(a;;x,t), 



where E(u;;a;,t) now denotes the monochromatic plane wave of Eq. (2). The integral over 
frequencies is easily evaluated for each of the regions of propagation, and the expression for 
the /? component of the electric field is 



Ep{x,t) 
Ep{x,t) 

Epix,t) 



V2 



(cos (3 + sin (3) , 



x<0 



COS p e CTc2 ^ no ' e ^ "o ' 



Eq 



+ sm pe "c^^ "e ' e = ^ "e 



< X < d 



Eq 

71 



cos/3e 



J^{x-ct-{l-no)df ^i'^{x-ct-{l-no)d) 



+ sin/3 e-^(^-^-('-'^«)'^)'eC(^-'=*-(i-"e)<i) 



d < X . 



(18) 



Even if this component is filtered out only after the signal has left the crystal, it is instructive 
to study its propagation in all the three regions, and we show a plot of the shape of this 
component for a particular choice of parameter values below. 

We first make a comparison between the peak position of the filtered gaussian signal and 
the corresponding position determined by the stationary phase argument. It is convenient to 
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change to dimensionless variables by measuring length in units of the crystal width d, and as 
coordinate in the direction of propagation we choose 

^=^ + (n-l), (19) 

so that outside the crystal the propagating waves are stationary in this coordinate. ^ = 
corresponds to the position midway between the gaussians of the ordinary and extraordinary 
waves after leaving the crystal, while the the incoming gaussian wave has the peak position 

^0 = n - 1 . (20) 

The peak position of the filtered signal, determined by the stationary phase argument, is 
under half waveplate conditions, 

= -^r + (n-l) 

2 cos/3 -sin/?' ^ ^ 

and the condition for superluminality, as discussed earlier, is then expressed as > ^o- 

The true peak position of the time evolved gaussian signal can now be determined by use 
of Eq. (18). For this signal we assume the central frequency u>o to satisfy the half waveplate 
condition 

^And= (2iV + l)7r. (22) 

We consider the squared relative amplitude / = (|i?^|/£'o)^ of the outgoing wave, which 
written as a function of ( has the form 

fiO = \ (cos^ exp(-/x2(C - \/^nf - sin/3 exp(-/i2(^ + \^nf^ . (23) 
The maxima (and minima) of /(^) are determined by the equation 

cos /3 (^ - ^An) e^'^^" - sin /? (C + ^ An) e'/^'^^" = , (24) 

with 

^=^, ^ = V^, (25) 

where £ measures the width of the incoming gaussian wave. 

We note that Eq. (24) reproduces the solution ^ = in the limit fi^^An 0. This 
means that the inequality yu^^An << 1 has to be satisfied if the position determined by the 
stationary phase argument should be a good approximation to the real peak position of the 
outgoing signal. If wc combine this condition with the condition for superluminality, > ^O) 
we have the constraints 

/ cos + sin B . , / — - — ; , 

^»\ n . ^ And>J2An{n-l)d. 26 

y cos p — smp 



8 



-1.5 -1 -0.5 



0.5 1 1.5 



Figure 4: Advanced and retarded signals. The figure shows the gaussian envelopes (the 
squared amplitudes /) of the fast (green) and slow (blue) components after leaving the crystal. 
For polarization angle (3 = Tr/4 there is destructive interference between the two components 
which gives rise to a symmetric signal with an advanced and a retarded maximum (red curve) . 

Thus, the width £ of the signal has to be comparable to or larger than the width d of the 
crystal in order for the time shift determined by the stationary phase argument to give a 
good approximation to the true time shift of the peak of the signal. Eq. (26) also shows that 
the closer P is to the value it/4, the larger the width i has to be. We further note from (23) 
that the amplitude of the filtered signal decreases exponentially with and therefore with 
the forward time shift r of the signal. 

In reality /(^) has two peaks rather than one, as follows from (23), since the function is 
defined by the partial cancellation of two shifted gaussians of opposite signs. This cancellation 
produces one advanced and one retarded peak. In Fig. 4 this is shown for /3 = 7r/4. For a 
wave with sharply defined value of k, this value of /? corresponds to the situation where 
the P component of the field vanishes. However, for a wave with a finite width in /c-space 
there is only a partial cancellation between the fast (ordinary) and the slow (extraordinary) 
components. The resulting wave is symmetric with respect to the advanced and the retarded 
peaks, as shown in the figure. For values of P smaller than it/ 4, there are still two peaks, but 
now there is an asymmetry, with the advanced peak dominating the retarded one. 

In Fig. 5 this situation is illustrated by plots of the function /(^) for the value /? = 0.2l7r 
and for three different values of the width i of the gaussian, one larger and two smaller than 
d. For the signal with smallest width the advanced and the retarded peaks are both clearly 
seen, however the amplitude of the retarded peak rapidly decreases with increasing width. 
One also notes that the forward shift of the advanced peak is close to the value for the 
broad signal with i > d, but is reduced below this value for the gaussians of smaller widths. 

Finally, in Fig. 6 we show the time evolution of the P component of the electric field, as 
the signal is transmitted through the crystal, in a case where the gaussian is sufficiently wide 
for the inequalities (26) to be satisfied. The width of the gaussian in this case is, £ = 4d, 
and value of the polarization angle is /3 = 0.2l7r. To enhance the time shift r in this plot, 
the degree of birefringence is increased to An = 0.5. The plot shows how the P component 
is suppressed under the propagation in the crystal, so that only a small component survives 
and appears as the outgoing signal. The position of the peak of this signal can be seen to 
be advanced relative to a gaussian signal that propagates with the speed of light through the 
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Figure 5: Outgoing signals for gaussian input signals of different widths. The three curves 
show the squared amplitude / as function of the dimensionless x-coordinate ^ for polarization 
angle (3 = 0.2l7r. For the (blue) curve with highest peak the width H of the gaussian is given 
by d/l = 2.6, with d as the width of the crystal. The corresponding values are d/i = 1.6 for 
the (green) curve of medium height and d/l = 0.6 for the broadest (red) curve. The position 
^0 indicates where the peak position would be if the gaussian signal moved with the speed of 
light, and if the signal moved with the speed determined by the stationary phase argument. 
These positions depend on the indexes of refraction which have here been chosen as n = 1.30 
and An = 0.15. 



crystal. 

Conclusion 

We have in this paper examined the phenomenon of superluminal propagation of electromag- 
netic waves through a birefringent crystal, an effect earlier discussed in Ref. [10]. The effect 
can be viewed as a two step process, where at the first step the crystal modifies the polar- 
ization of the propagating wave, due to different velocities of the ordinary and extraordinary 
components, but since there is no damping or reflection, the signal travels with subluminal 
velocity through the crystal. However, at the next step, when the signal is sent through a 
polarization filter, the effect of the filter is to suppress the signal in such a way that the peak 
of the outgoing signal will appear as being shifted forward in time. 

The amplitude of the filtered signal has been examined here as function of the frequency 
and polarization angle for an incoming monochromatic plane wave. Close to half waveplate 
condition, a stationary phase analysis indicates that the filtered signal may be substantially 
shifted forward in time, in agreement with a similar analysis in Ref. [10]. The implication for 
a modulated wave is that, for certain parameter values, the position of the peak of the filtered 
component will be shifted to an advanced position relative to a freely propagating wave. This 
is what in [10] is interpreted as a superluminal transit of the signal through the crystal. 

To examine this effect further we have considered incoming waves with gaussian envelopes 
of varying widths. The effect of the filter on the transmitted wave will in general be to reduce 
the signal in such a way that a two-peaked signal will appear, with an advanced and a retarded 
part. For parameter values corresponding to superluminal transit time the advanced peak will 
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Figure 6: Propagation of the /3-component of the gaussian wave through the birefringent 
crystal. The figure shows the envelope of the component at four subsequent times, ti, t2, 
and t^. The polarization angle is /3 = 0.2l7r and the width is i = Ad. The region where the 
crystal is located is indicated by the shaded (blue) area. The propagating wave of polarization 
l3 is shown by the red curve, while the green curve is included for comparison to show the 
propagation of the same signal in empty space. The peak position of the transmitted signal, 
which is indicated by the dotted vertical line, clearly shows that the peak has an advanced 
position relative to the peak of the signal in empty space. In this plot the following values 
have been chosen for the indexes of refraction, h = 1.35 and An = 0.5. 

be the dominating one. However, a good agreement between the stationary phase analysis of 
the peak position and the true peak position of the gaussian signal depends on a sufficiently 
broad gaussian envelope in the direction of propagation. We find that if the width of the 
signal is small compared to the width of the crystal, the peak will be located behind the 
position determined from stationary phase analysis and the filtered wave in this case will 
also show both the advanced and retarded peaks. If instead the width of the gaussian is 
much larger than the crystal width, the peak position will be close to that determined from 
stationary phase and it may be substantially advanced relative to signal moving with the 
speed of light. However, the broadness of the signal means that it is not strongly located 
around this position. 

The stationary phase analysis shows that the situation discussed here is similar to other 
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cases, where reduced transit time can be seen as a consequence of the anomalous dispersion 
of the transmitted signal. A particular such case is quantum tunnehng, where the outgoing 
wave may appear with reduced amphtude but advanced position [15, 16, 17]. However, the 
present case is, in a sense, more transparent than the other cases due to the separation of the 
effect into two parts. Thus, at the first step, when the wave propagates through the crystal, 
there is no faster-than-light propagation, but only a change of polarization. At the next step, 
the polarization filter produces the advanced signal by separating out a special superposition 
of the two subluminal waves. That means that the superluminal effect, in the present case, 
can be viewed as due to fast transit of interference patterns through the crystal. Viewed in 
this way the effect seems close to other cases where velocities larger than c can be obtained 
by interference. A particular case is the superluminal X-wave in free Maxwell theory, which 
is created by interference between monochromatic plane waves [18, 19]. 
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